Momentum dependence in the dynamically assisted Sauter-Schwinger effect 
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Recently it has been found that the superposition of a strong and slow electric field with a weaker 
and faster pulse can significantly enhance the probability for non-perturbative electron-positron 
pair creation out of the vacuum - the dynamically assisted Sauter-Schwinger effect. Via the WKB 
method, we estimate the momentum dependence of the pair creation probability and compare it 
to existing numerical results. Besides the theoretical interest, a better understanding of this pair 
creation mechanism should be helpful for the planned experiments aiming at its detection. 



PACS numbers: 12.20.-m, ll.15.Tk. 
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I. INTRODUCTION 

Quantum electrodynamics (QED), the theory of 
charged particles such as electrons and positrons inter- 
acting with the electromagnetic field, offers many inter- 
esting phenomena beyond standard perturbation theory 
(i.e., expansion into powers of the fine-structure constant 
aQED). One prominent example is the Sauter-Schwinger 
effect describing the ereation of electron-positron 

pairs out of the QED vacuum by an electric field via tun- 
neling. For a constant electric field E, the lowest-order 
pair creation probability scales as 
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where m is the mass and itg the charge of the elec- 
trons and positrons. Since this expression does not 
permit a Taylor expansion in this phenomenon is a 
non-perturbative effect. Because the characteristic field 
strength Es (which is also called the critical field) is very 
large ^ 10^®V/m, this fundamental QED prediction has 
not been conclusively experimentally verified yet - in con- 
trast to electron-positron pair ereation in the perturba- 
tive multi-photon regime, see, e.g., Q. 

Furthermore, in spite of many efforts, see, e.g., (5l-[l0|. 
our understanding of this non-perturbative effect for non- 
constant electric fields E is still far from complete. For 
example, recently it has been found [Tl| that the superpo- 
sition of a strong and slow electric field with a weaker and 
faster pulse can significantly enhance the pair ereation 
probability, see also [l^ - IT^ . In the following, we are go- 
ing to present an analytic estimate of the momentum de- 
pendence of this enhancement mechanism via the WKB 
approximation (for more details, see also [l^). These 
findings should be relevant for the planned experiments 
with optical lasers or XFEL (or a combination of both) , 
see, e.g., [l^, which might facilitate the first conclusive 
observation of this effect. 



II. BASIC FORMALISM 

We consider quantized electrons and positrons in an 
external (classical) electric field E with no magnetic field 
present B = 0. The field E is supposed to be purely time- 
dependent E = E{t)ex and sub-critical E <^ Es- Thus 
we can neglect the spin of the electrons and positrons and 
use the Klein-Foek-Gordon equation in temporal gage 
E = A where A is the vector potential {h = c = 1) 



[V - iqA{t)Y 



$ = 0. 



(2) 



After a spatial Fourier transform, we get [1] 



0. 



(3) 



We see that the momentum fej_ transversal to the electric 
field can be absorbed by a re-definition of the effective 
mass rn? — >■ m^g = k\ ~\- m? . Thus, we omit it in the 
following formula for brevity. 

The above equation is formally equivalent to a har- 
monic oscillator with a time-dependent potential 



(4) 



If we replace t by x, <j)[t) by ili[x), and ^l'^[t) by 
2m[E — V{x)], this equation has the same form as a one- 
dimensional Sehrodinger scattering problem with energy 
E and potential V{x). Therefore, a solution which ini- 
tially behaves as (p^^it) = exp{— iwfct} will finally evolve 
into a mixture of positive and negative frequencies 

= "fc exp{-iujkt} + Pk exp{+iwfci} , (5) 

where the Bogoliubov coefficients and jSk are related 
to the reflection R and transmission T amplitudes in the 
one-dimensional scattering theory picture via otf, = 1/T 
and pf, = R/T. The probability for electron-positron 



pair creation is given by 



III. RICCATI EQUATION 
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For slowly varying and sub-critical electric fields E(t), 
the Bogoliubov coefficients oif, and /3fc can be derived via 
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the WKB approximation. To this end, let us cast Eq. ^ 
into a first-order form 



d 
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M v. (6) 



Since usual scalar product is not conserved by this evolu- 
tion equation, it is useful to introduce the inner product 



{v\v') = ~~i{vlv'2 - V2V[) 



(7) 



which is just the Wronskian of the original Eq. Q. Thus, 
the inner product of two solutions v and v' is conserved 



dt 



{v\v') = 0. 
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As the next step, we expand the solution v{t) of Eq. ^ 
v{t) = a(t)e*'^(*)t;+(t) + /3(t)e-*'^(*'t;_(t) , (9) 
into instantaneous eigenvectors v±{t) of the matrix 

M ■ v±{t) = ±in{t)v±{t) . (10) 



With the normalization v± = (1, ±if2)-^/A/2J7, we find 
{v+\vjf-) = 1, {v-\v^) = —1, and («+|«_) = 0. Finally, 
a{t) and /3(t) are the instantaneous Bogoliubov coeffi- 
cients, where we have separated out the WKB phase 



ip{t) = / dt'Q{t') . 



(11) 



This uniquely determines the evolution of a{t) and /3{t) 
which we can obtain by inserting the expansion ^ into 
the equation of motion (jS]) and projecting it with the 
inner product ([7|) onto the eigenvectors v±(t) 

where we have used {v+\v^) = = 0, as well as 

{v^\v+) = n/{2n) and = -i7/(2^]). In terms 

of the reflection coefficient R{t) = l3{t)/a{t), we get 

m = 1^ - i?^(Oe-^'^(*)) , (13) 

which is known as Riccati equation, see also 0, [3| • 



IV. WKB METHOD 



The Riccati equation ([13)) is still exact but unfortu- 
nately non-linear. The WKB approximation is based on 
the assumption that the rate of change of n{t) is much 
slower than the internal frequency n{t) itself. In our case. 



n^(t) ^ [k^ - qA{t)]^ + m'^ 



(14) 



this is satisfied if the strength E and the rate of change uj 
of the electric field E{t) are small compared to the mass, 
i.e., for E <^ Es and u) ^ m. In this limit, the phase 
factors e^^"^ are rapidly oscillating and the magnitude 
of R can be estimated by analytic continuation to the 
complex plane. In terms of the phase variable Lp, the 
Riccati equation (|13p reads 



dR{Lp) 
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(15) 



Now analytic continuation to the upper complex half- 
plane (f — > ^(p + i'^ip = 5R(/3 + ix shows that R becomes 
exponentially suppressed R ^ e~^^. The degree of this 
suppression depends on the point where the analytic con- 
tinuation breaks down. Since e^^"^ is analytic every- 
where, this will be determined by the term In $7. Typi- 
cally, one can go into the upper complex half-plane until 
one hits the first zero of il at i*, i.e., fi(i*) = 0. These 
points in the complex plane are analogous to the clas- 
sical turning points in WKB. Consequently, we find 



i?(i t oo) = i? - e 



exp{-25[(^(i,)]} . (16) 



So far, this is only an order-of-magnitude estimate. How- 
ever, it can be shown that this expression becomes 
exact (under appropriate conditions) in the adiabatic 
limit (roughly speaking, m f oo), i.e., that the pre- factor 
in front of the exponential tends to one. 

In case of more than one turning point, the one with 
the smallest x*j i-e., closest to real axis (in the complex 
(y9-plane) dominates. For multiple turning points with 
similar x* i there can be interference effects [3, 0] • 



V. DOUBLE SAUTER PULSE 

Now we arc in the position to apply the above method 
to the double Sauter pulse studied in [ll| 
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Note that the Dirac equation in the presence of a single 
pulse (e.g., E2 = 0) of that type can be solved exactly. 
This solution has already been used by Sauter [l[ even 
though with t and x interchanged. 

The first term on the r.h.s. describes a strong and slow 
electric field profile while the second term corresponds to 
a much weaker and faster pulse with wi W2 w and 
E2 El <^ Es- The Keldysh ^ parameters of the two 
pulses are supposed to be small 71 = mtoi/ (qEi) <^ 1 
and large 72 = muj2/{qE2) ^ 1, respectively, while the 
combined Keldysh parameter introduced in [ll| 



mLij2 



(18) 



is of order one. For time scales where the second (fast) 
pulse contributes, we may approximate the first (slow) 
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pulse by a constant field Ei due to wi ^ a;2. Thus the 
vector potential becomes 



A{t) Ri Eit + — tanh(a;2i) , 

UJ2 



(19) 



and the condition fi(t*) ~ for the turning points reads 

1^2^* + tanh(w2i*) = 7c ( « + — ) ■ (20) 
El \ mj 

Due to the periodicity of the function tanh(a;2i)i there 
are infinitely many solutions. However, we are mainly 
interested in those close to the real axis - which have the 
smallest value of x* ■ For Ei ^ E2, we may approximate 
the two most relevant solutions t* and t' via 
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Note that this approximation breaks down when 7c ap- 
proaches 7r/2, so wc will assume 7c > 7r/2 in the follow- 
ing, see also [ll|. The first solution t* is basically the 
same as for the slow pulse alone (i.e., E2 =0), while the 
second one t'^ is obviously tied to the fast pulse. 

Again using E^ >• E2, we may derive the associated 
exponents x* ^-i^d x'*- For the normal solution t*, we 
reproduce the usual value given solely by the slow pulse 
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With Pe+e- ~ l^^l - e-4x., we verify Eq. ^. 

This result is independent of because we have approx- 
imated the slow pulse by a constant field. The anomalous 
solution t'^, on the other hand, yields 



^5 / — 
2qEi V \2j. 
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with f{x) = x\/l + + arcsinh(a;). In this form, the 
dependence on kx is perhaps not so obvious. Thus, let 
us Taylor expand this expression in powers of kx 
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The momentum- independent result in the first line was 
already derived in [Tlj via the word-line instanton 
method, but this method did not yield any information 
about the momentum dependence in the second line. 



VI. DISCUSSION 

For the dynamically assisted Sauter-Schwinger ef- 
fect [ll| . we estimated the momentum dependence in 
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FIG. 1: (Color online) Comparison of analytical and numeri- 
cal results for the (differential) pair creation probability |_Rp 
depending on the momentum p = k^c/meS- The solid (blue) 
curve shows the numerical data from (l3 | (kindly provided 
by the authors) for Ei = Ss/4, cji = 10*eV, E2 = -Bi/10, 
and (jj2 = 5 X 10'''eV. For the same values, the function 
j-Rp — |c, exp{— 2(^(f*)}-|-c', exp{— 2(/3(f'^)}p representing the 
WKB result is plotted (dashed purple curve), where and 
t', are the (exact) WKB turning points and ip is the complex 
phase in Eq. The pre- factors are chosen as c, = 1.4 and 

= 0.2. Finally, the dotted (brown) curve describes the case 
-B2 = of a single pulse (i.e., the background). 



Eqs. (221) and ^ via the WKB method. The strong 
and slow pulse alone is represented by the normal turn- 
ing point with x* i'^ Fq. (j22p . It has a broad mo- 
mentum distribution {Akx)'^/m'^ ^ l/7i ^ 1 which can 
be explained by the uncertainty At ^ l/wi of the cre- 
ation time: Particles that are created earlier have more 
time to be accelerated by the st rong electric field than 
others which arc produced later [13|. The dynamically 
assisted effect, on the other hand, corresponds to the 
anomalous turning point t'^ determined by the weak and 
fast pulse. Thus, the creation time is far less uncertain 
At ^ 1/1^2 and - as one would expect - it has a much 
narrower width, even though it is probably hard to guess 
the precise scaliiig {AkxY ~ qE^^^"^ — tt^ beforehand. 
Consistent with [ll|, this window closes for 7c J, 7r/2, af- 
ter which the enhancement disappears. For 7c < 7r/2, 
the normal turning point t» dominates x* < x'* • 

Let us compare our findings with the numerical results 
in [1^ obtained via the quantum kinetic approach. In 
Fig. [TJ we plot the numerical data from Fig. 5a of [l^l 
as well as \R\^ = \c* exp{-2ip{t^)} + c'^ exp{-2(^(<'J}P, 
where ip is the complex phase in Eq. ([TT]). The two WKB 
turning points t^ and t'^ can be obtained exactly, but our 
WKB method does not yield the pre-factors c, and c'^, 
which can also depend on p in general. Neglecting this 
p-depcndence, we may estimate the normal pre-factor c* 
by comparison to the case of a single pulse (i.e., E2 = 0) 
which facilitates an analytic solution. The anomalous 
pre-factor c'^,, on the other hand, was chosen (fitted) to 
match the numerical results. After that, the agreement 
between the analytic and the numerical results is surpris- 
ingly good ~ given that the employed values Ei — Es/i, 
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wi = lO'^eV, E2 = Ei/IQ, and a;2 = 5 x lO'^eV do not 
satisfy our underlying assumptions E2 Ei <^ Eg and 
uji ^ UJ2 m very well. Note that the differenee be- 
tween c* and c'^ also indicates that we are not deep in 
the adiabatic limit (m f 00). Furthermore, one should be 
very careful with the order of the various limits in this 
multiple-scale problem. For example, the adiabatic limit 
(m t 00) does not commute with the limit E2/E1 | 
since the dynamically assisted Sauter-Schwinger effect 
given by c'^ should vanish for E2 = 0- Finally, the os- 
cillations visible in Fig. [1] (and Fig. 5a o f [l4| ) can be 
explained nicely by interference effects 0, Qof the two 
turning points and t'^. The interferences are most 
pronounced where the two contributions c* exp{— 2(p(t*)} 
and c'^, exp{— 2i^(t'^)} are equally strong, which happens 
around ~ irrics in this case. 

VII. OUTLOOK 

It might be interesting to generalize the above findings 
to other pulse profiles such as 

E{t) = Eihiuit) + E2h{^2t) , (25) 

with /i(0) = h{^) = 1 as well as wi ^ W2 ^ w and 
E2 <^ El <^ Es- For a broad class of functions, for 



example f2{x) = 1/(1 -|- x^)^, we expect to obtain qual- 
itatively the same picture as discussed above. Again, 
the normal turning point will have basically the same 
value as in Eq. (j2ip and the anomalous turning point t'^, 
will be very close to the singularity of /2 , in our example 
t'^ « i/uj2- As before, for small values of the combined 
Keldysh parameter 7c in Eq. (|18|) . the contribution of the 
normal turning point dominates - but if 7c exceeds a 
critical value of order one, the anomalous turning point 
t'^ becomes stronger than the normal one and we get dy- 
namically assisted pair creation. 

For other profiles, such as f2ix) = exp{— a;^}, however, 
the anomalous turning point t'^ logarithmically depends 
on the ratio Ei /E2 and thus the mechanism of dynami- 
cally assisted pair creation (including the threshold value 
for 7c) will also depend on E1/E2, in contrast to the case 
considered above. 
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